Theoretical Evidence for the Berry-Phase Mechanism of Anomalous Hall Transport: 

First-principles Studies on CuCr2Se4_^.Br2^ 
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To justify the origin of anomalous Hall effect (AHE), it is highly desirable to have the system 
parameters tuned continuously. By quantitative calculations, we show that the doping dependent 
sign reversal in CuCr2Se4-a;Bra:, observed but not understood, is nothing but direct evidence for 
the Berry-Phase mechanism of AHE. The systematic calculations well explain the experiment data 
for the whole doping range where the impurity scattering rates is changed by several orders with Br 
substitution. Further sign change is also predicted, which may be tested by future experiments. 

PACS numbers: 75.47.-m, 71.20.-b, 72.15.Eb 



In spite of the wide applications of anomalous Hall ef- 
fect (AHE) to characterize ferromagnetism, its origin has 
been a controversial subject since its discovery more than 
a century ago [1]. The k-space gauge fields, known as 
the Berry curvature, exist ubiquitously in Bloch bands 
where time reversal symmetry is broken, giving rise to 
an intrinsic anomalous Hall effect (AHE) in ferromag- 
nets 0. This intrinsic effect was originally derived by 
Karplus-Luttinger fifty years ago based on a linear re- 
sponse theory [3], but was disputed ever since and until 
recently, extrinsic mechanisms of skew scattering and side 
jump were usually invoked Inspired by the new un- 
derstanding from the Berry phase connection @, @, 0]j 
a number of quantitative studies have been successfully 
carried out in recent years P, loj , finding that the Berry- 
phase mechanism is important in various materials. How- 
ever, theoretical understanding of the condition for such 
importance is far from clear, despite a large number of 
theoretical analysis based on model Hamiltonians [lo[. 
To fully explore the importance of the Berry-phase mech- 
anism, it is highly desirable to have a systematic study 
of real materials in comparison to experiments when the 
system parameters are tuned continuously. 

In this paper, we report systematical first-principles 
calculations on doping-dependence of the intrinsic AHE. 
Our material of choice is the ferromagnetic spinel, 
CuCr2Se4, one of the parent compounds of a wide class of 
colossal magnetoresistive chalcospinels. It is well known 
for its high Curie temperature (Tj. = 450 K) and large 
room-temperature magneto-optic Kerr effect, with great 
potential for spintronics applications The experi- 

mental measurement of AHE in this compound has been 
recently carried out by Lee et al. [l^ , where they are able 
to tune the scattering rate by 70 folds with Br substitu- 
tion of Se. Our quantitative calculations well explain the 



experimental AHE data over the whole doping range with 
reasonable accuracy based on the Berry-Phase mecha- 
nism. In particular, we reveal that the sharp sign change 
in the doping dependent anomalous Hall conductivity, 
which was observed in the experiment but not discussed 
explicitly, is a direct evidence for the Berry-Phase mech- 
anism of AHE. The sign change is due to a large patch 
of high Berry-curvature in the band structure. In addi- 
tion to explaining this experiment, our calculations also 
extend to the case of hole doping, urging further experi- 
ments on the spinel system. 

The spin-polarized ground state of CuCr2Se4 has been 
calculated by the Xa method Il3| and by the lin- 
earized muffin-tin orbital method [ill , [l^ . In this work, 
the relativistic electronic structure is calculated self- 
consistently using the full-potential linearized augmented 
plane- wave method with generalized gradient approxima- 
tion (GGA) [l5[ for the exchange-correlation potential. 
We use the experimental lattice constant, and the mufhn- 
tin radius Rmt = 2.1,2.4,2.3 Bohr for Cu, Cr, and Se 
atoms, respectively. The convergence of present calcula- 
tions has been well checked. 

Figure [1] shows the calculated total and projected den- 
sity of states of the parent compound CuCr2Se4, where 
Cr atoms occupy the octahedral sites, and Cu atoms 
occupy the tetrahedral sites. To understand the com- 
plicated electronic structure, we consider the compound 
as a combination of two parts: the tetrahedral (CuSe4) 
clusters (in the 6-f nominal valence) which are arranged 
periodically in the crystal space of the diamond struc- 
ture, and the Cr atoms (with 3+ nominal valence) in the 
interstitial sites of the CuSe4 (diamond) crystal frame- 
work. As shown in Fig.[Tl the electronic states of (CuSe4) 
is almost non-spin-polarized (the slight polarization will 
be discussed later). The Cu is nearly in the Cu+ va- 
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FIG. 1: The calculated total and projected density of states 
of CuCr2Se4. The Fermi level locates at energy zero. 



lence state, whose 3d orbitals are almost fully occupied 
and are away from the Fermi level. The electronic states 
around the Fermi level mostly come from the Cr 3d and 
Se 4p states. The 3d states of Cr^+ are exchange split by 
about 3.0 eV, giving rise to the high spin configuration 
itlle"H^^e"^) with 3.0 /is/Cr local moment. Here the Cr 



3ci— Se 4p hybridization is an essential factor to form the 
final electronic structure. First, the hybridization will in- 
duce holes in the Cr states, resulting in reduced local 
moment and enhanced valence (Cr^+''). This is evident 
from the slightly non-occupation of Cr t^^ states around 
the Fermi level (see Fig. [T]). Second, the hybridization 
leads to the negative spin polarization of itinerant Se Ap 
states (anti-parallel to the spin moment of Cr), which is 
crucial for the AHE in this compound. Finally, the hy- 
bridization stabilizes the ferromagnetic ground state and 
contributes to the high Curie temperature as discussed 
for SrFeMoOe and (GaMn)As [i3|. The calculated total 
moment is 5.1 /is/f.u. for the parent compound, which 
is in good agreement with the experimental value of 5.2 

The Br substitution will introduce additional electrons 
(in addition to increasing disorder) due to the reduced 
negative valence of Br compared with Se. It is justi- 
fied by the following facts that the electronic structure 
with doping can be described by the rigid-band shift {i.e. 
changing doping is equivalent to sweeping the Fermi en- 
ergy) without losing the main physics for our purpose. 
(1) It was reported IJ] that the Br substitution only 
affects the Curie temperature, but does not affect the 
ferromagnetic ground state dramatically. (2) By 25% 
substitution (a; = 1.0), the lattice parameter changes 
only by 0.7% [l3|- (3) As a self-consistent check, the ob- 



TABLE I: The calculated spin, orbital and total moments of 
CuCr2Se4_2!Br3: in unit of fiB. 
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tained electronic structures with the rigid-band approx- 
imation is used to calculate the magnetic moments and 
gives results in good agreement with experimental data. 
As shown in TablelH the calculated total moment per f.u. 
increases monotonically from 5.1 hb for x = 0.0 to 6.0 iib 
for X — 1.0, while the experiment shows an increase from 
5.2 fiB to 6.0 fiB The calculated spin and orbital 

moments of each atom also agree well with the results by 
X-ray magnetic circular dichroism studies [lH. The or- 
bital moment of Se sites mainly comes from its 4p states 
due to the spin-orbit coupling. 

The intrinsic anomalous Hall conductivity can be eval- 
uated from the linear response theory using the standard 
Kubo formula [9]] 
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where \ipnk) is the eigenstate with eigenvalue Enk = 
hujnk, Vx and Vy are the velocity operators, fnk is the 
Fermi-Dirac distribution function, and <5 is a small pa- 
rameter representing the finite life-time broadening of the 
eigenstates. The $!„(&) is a vector in the fc-space, and 
can be related to the Berry curvature of the Bloch state 
in the clean limit {5 = 0), i.e., r2„(fc) = Ini(VfcW„fc| x 
|VfcU„fc) with Unk being the periodic part of the Bloch 
wave function. 

Figure [2] shows the calculated intrinsic 

(7xy cL func- 
tion of the doping x. Let us first consider the calculated 
curve in the clean limit {5 ~ 0, open circles). It is obvi- 
ous that Uxy is highly non-monotonic and changes its sign 
twice between x = 0.0 and 0.5; it starts with a positive 
value at a; = 0, then changes its sign to negative around 
X — 0.1, and again to positive for x > 0.3. The places 
where axy changes sign with varying x, although appear 
arbitrary, are in good agreement with the experimental 
data (the square-cross points in Fig. (2]). While a quan- 
titative comparison with the experiment of the overall 
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FIG. 2: (Color online). The anomalous Hall conductivity a^y 
as function of doping x in CuCr2Se4-a;Bra; . For the theoretical 
results, the open circles are a^y for the clean limit (5=0), 
and the triangles are a^y with finite 5 in the Kubo formula, 
where the doping dependent S is determined from the ab-mitio 
calculated plasma frequency (shown in the lower panel) and 
the experimental longitudinal resistivity (see the text part 
for details). The square-crosses are experimental results from 
Ref. fii. 



behavior of a^y needs further analysis (as addressed be- 
low), such an agreement is a striking result, considering 
the fact that the calculations were done systematically 
for the whole doping region without adjustable parame- 
ters. In the experimental analysis [l^ . the sign of (Jxy 
is dropped; only its absolute value is taken into account. 
Our result here, however, shows that the sign of Uxy is 
important, and the sign change of a^y with varying dop- 
ing x is a natural result of the Berry-phase mechanism 
of AHE. 

To make the quantitative comparison, we need consider 
the effect of the finite life-time of the eigenstates. The 
simplest way to do this is to assume the diagonal form 
of electron self-energy and to use a single parameter 5 
instead (in the Kubo formula), thereby neglecting the 
vertex correction due to impurity scattering. It is worth 
to note that the doping dependent 5 in our approach is 
not adjustable parameter but instead it is determined 
from the relaxation time t = h/5 ~ l/{eoUjpp), where 
p is the longitudinal resistivity, adopted from Ref. [l^ . 
and the plasma frequency uip is calculated from the band 
structure by 
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The plasma frequency is actually the measurement of the 
ratio between the number of band carriers n* and the 
effective mass of electrons m*, according to the relation 
ujp = n*e^/(£om*). The triangle points in Fig. [2] are the 
theoretical values of dxy after considering the effect of 
relaxation. It is now obvious that the calculated intrinsic 
axy is in quantitative agreement with experimental data, 
especially in the region around the sign change {x = 0.3). 
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FIG. 3: The calculated band structure of CuCr2Se4_a;Bra; 
(upper panel) and the sum of Berry curvature over the oc- 
cupied bands 51^ (fe) = X^„/nfcn^(fe) for two characteristic 
Fermi level positions corresponding to doping x — 0.21 and 
0.37 respectively. The lower panels show the Q^{k) in a spe- 
cial plane of the BZ for different doping x. 



The doping-dependent sign change of axy was observed 
before in other ferromagnetic alloys, such as Nix-xF^x, 
Au-Fe and Au-Ni alloys [3]. The phenomenological the- 
ory [l^ relates the sign change to the change of effec- 
tive spin-orbit coupling with varying chemical potential. 
Here our numerical calculation indicates that the sign 
change in CuCr2Se4_a;Bra; is attributed microscopically 
to the topological nature of electronic bands in the Berry 
phase picture. From simplified two-band mode, it is un- 
derstood that the sum of Berry curvatures over the oc- 
cupied bands n^(fc) = J^nf^k^nif^) is spiky and the 
sign change occurs near the degenerate or band crossing 
points, which act as magnetic monopoles in the momen- 
tum space 0. As a result, by summing over the Brillouin 
zone (BZ), axy is typically a non-monotonic function of 
chemical potential, and exhibits sharp fluctuation. This 
is the case for CuCr2Se4„2:Br2; as shown in Fig. [3l A 
similar behavior was also observed in the 2-dimensional 
(2D) systems, such as the sign change in the quantum 
well structure 
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However, we note that the higher 
dimensionality in the present system makes the problem 
quite different. In the 3D case, the single band crossing 
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point cannot contribute enough weight to the sign change 
of axy due to the 3D (instead of 2D) integration of BZ. 
To get enough weight, a high density of states near band 
crossing points (or near degenerate points) are necessary 
(for example, the insert in Fig. [3]). Due to the presence 
of band dispersion, it is generally hard to have all those 
band crossing points occupied (or unoccupied) at each 
fixed chemical potential, which leads to lower possibility 
for the sign change of a^y in 3D than in the 2D case. 
On the other hand, CuCr2Se4_2,Brj, is an isotropic 3D 
system where sharp sign changes of axy are observed. 
Actually, the sign changes in CuCr2Se4_2,Br2, are neither 
from simple band crossing nor from the high symmetric 
points of the BZ. As shown in Fig.[3l the dominant nega- 
tive Berry curvature for x — 0.21 (the valley of Uxy) and 
the positive Berry curvature for x = 0.37 (the peak of 
axy) are located at different regions of the BZ. We have 
tried to use an effective Luttinger Hamiltonian (fitted 
from our electronic structure calculations) to study the 
system, but the sign changes cannot be reproduced even 
qualitatively. This indicates that in realistic materials 
the accurate first-principles calculations are important. 

In conclusion, the doping-dependent AHE in 
CuCr2Se4_a;Bra; is investigated by ab initio calcu- 
lations, and analyzed according to the Berry phase 
picture. The good agreement between experimental 
and numerical results provide strong evidence for the 
Berry-phase mechanism of AHE, even when the impu- 
rity scattering rates is changed by several orders. The 
disorder (extrinsic) contributions, which may also be 
related to the non-zero Berry curvature (20l | , can change 
the magnitude of our calculated AHE quantitatively, 
but they are not expected to affect such features as the 
sign change qualitatively. To further verify our results, 
we point out the following two aspects which can be 
justified experimentally. (1) Additional sign change is 
predicted from our calculation. As shown in Fig. [51 
by negative doping (hole doping), we predict that axy 
changes its sign from positive to negative. The hole 
doping can be realized experimentally by doping As 
instead of Br. (2) The experimentally observed Nernst 
effect 2l| in the same compound CuCr2Se4_a;Br^ can 



be also checked from the present picture (22i] . 
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